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ABSTRACT 

We investigate the mass loss from a rotationally distorted envelope following the early, 
rapid in-spiral of a companion star inside a common envelope. For initially wide, 
massive binaries (M1+M2 = 20 M@, P ~ 10 yr), the primary has a convective envelope 
at the onset of mass transfer and is able to store much of the available orbital angular 
momentum in its expanded envelope. Three-dimensional SPH calculations show that 
mass loss is enhanced at mid-latitudes due to shock reflection from a torus-shaped 
outer envelope. Mass ejection in the equatorial plane is completely suppressed if the 
shock wave is too weak to penetrate the outer envelope in the equatorial direction 
(typically when the energy deposited in the star is less than about one-third of the 
binding energy of the envelope). We present a parameter study to show how the 
geometry of the ejecta depends on the angular momentum and the energy deposited 
in the envelope during a merging event. Applications to the nearly axisymmetric, but 
very non-spherical nebulae around SN 1987A and Sheridan 25 are discussed, as well 
as possible links to RY Scuti and the Small Magellanic Cloud object R4. 
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1 INTRODUCTION 

The common-envelope (CE) phase is one of the most impor- 
tant and least und erstood pha s es of stellar evolution. Orig- 
inally proposed by iPaczvriskil fll9T6T) to explain the origin 
of short-period binaries with compact objects, it can also 
significantly alter the evolution of systems in which the en- 
vel ope remains bound, leaving an atypical single star (see 
e.g. Ilvanova fc Podsiadlow ski 2003). For example, it is now 
widely believed that the unusual properties of the progeni- 
tor of SN 1987A are due to a binary merger some 20 ,000 yr 
before the explosion iPodsiadlowski fc Ivanovall2003t) . 

In this paper we are interested in the case where the 
primary initiates mass transfer either when crossing the 
Hertzsprung gap (so-called early case B mass transfer) or 
later as a red supergiant (late case B/C mass transfer). In 
the latter case, the primary has already developed a deep 
convective envelope and mass transfer is dynamically unsta- 
ble if the mass ratio exceeds a critical value, leading to a 
common-envelope and spiral-in phase. 

Early case B mass transfer initially occurs on the ther- 
mal timescale of the mass donor and is dynamically stable; 
but the secondary may not be able to accrete all of the trans- 
ferred mass, and this may also lead to a common-envelope 
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system and possibly the merging of t he system (e.g. 
ll994tlWeUstein. Langer fc Braunlbool . 

While the common envelope maintains co-rotation with 
the embedded binary, orbital angular momentum is effi- 
ciently transferred from the binary orbit to the envelope, 
where most of the initial orbital angular momentum is stored 



Lorb = 6.60 x 10 54 gcm 2 s _1 A 1 ^ Mi 5 M 5 M~ 



1/2 



(1) 



where A2500 is the orbital separation in units of 2500 Rq , 
M15 and M5 are the masses of the primary and the sec- 
ondary in units of 15 Mq and 5 Mq (as indicated by the sub- 
scripts) , respectively, and M20 = Mi + M2 is the total mass 
in units of 20 Mq . This phase may last perhaps for a few 
decades and ends when the envelope becomes differentially 
rotating. The subsequent rapid plunge-in of the secondary 
then drives significant envelope expansion and the ejection 
of at le ast some of the envelope llMe^ej^^J^ev^r^Iofrneister] 
1 19791 ISandcmist et all Il998t iTaam fc Sandauisd l200of 
Podsiadlowski 2001; Ivanova & Podsiadlowski 2003). If the 
envelope is not completely ejected in this phase (the case of 
interest in this study), the spiral-in continues and now be- 
comes self-regulated where all the energy dissipated by the 
further orb ital decay is transported to the surface and radi- 
ated away jMever fc Mever-Hofmeisterlll979l : iPodsiadlowskj 
2001). A second phase of mass loss may result from a nu- 
clear flash that may occur du ring the final core merger 
ilvanova fc Podsiadlowskill2003l) . 



© 2005 RAS 



2 T. Morris and Ph. Podsiadlowski 




Figure 1. Particle plots showing specific angular momentum (increasing outwards) and angular velocity for the three cases of angular 
momentum L = 0.235, 0.588, 0.817V GM 3 R from left-hand panel to right-hand panel. The solid curves give the critical specific angular 
momentum j c while the dashed curves show the solid body profiles j = uir 2 where u = 1. All quantities are given in 

code units. 




Figure 2. Distribution of SPH particles in the meridional plane immediately after the spin-up for the three cases of angular momentum 
L = 0.235, 0.588, 0.817VGM 3 R from left-hand panel to right-hand panel. For other properties see Tabled] 



Previous one-di mensional numerical simulations by 
IPodsiadlowskil 1)200 j) have shown that significant mass loss 
may occur even when most of the en velope remains bound . 
In t he three dimens i onal m odels of iLivio fc Sokerl dl98ct) 
and ISandauist et alJ (1199a) . most mass loss occurs in the 
orbital plane of the binary. However, these authors consid- 
ered the case where most/all of the envelope was ejected. In 
this study, we consider the less energetic case appropriate 
for a merger. As we will show, in this case mass loss may 
preferentially occur at mid-latitudes and be suppressed in 
the equatorial direction if the energy deposited is less than 
about one-third of the binding energy of the envelope. In sec- 
tion 2 we outline our numerical method and in section 3 we 
present the main results of our study and their dependence 
on the input parameters. In section 4 we apply these results 
to observed systems, in particular SN 1987A and Sheridan 
25. 



2 NUMERICAL METHOD 

We model the common envelope as a condensed polytrope 
with adiabatic index 7 = 5/3 with a central point mass 
which contains two-fifths of the system mass . Assuming 



spherical symmetry and hydrostatic equilibrium initially, we 
obtain the radial density profile by integrating the dimen- 
sionless equations, 



dp _ _ .m _ 2 - 7 = f G Af c 2 orc 7 
dz z* P \lKR£*? 

dm 2 - 

— — = 47TZ p 

dz 

with the inner boundary conditions 

m(z = 1) = 1 

and 

p(z = 1) = (yA)^ I S(d ■ 



-39 



8tt 



(2) 
(3) 

(4) 
(5) 



where z 
The free 
boundary 
where Z 
presented 
ISpringeT~ 



ity and g 
dynamics 



= r/Rcov C , m = M(r)/M corc and p = pM^/Rl. 
parameters A, 8 are determined from the surface 
conditions p[z = Z) = 0, fh(z — Z) — M*/M cor o 
= R t /Rcoic- For all hydrodynamical simulations 
in this paper, we us e the GADGET code of 
Yoshida fc White! teOOll) . which implements grav- 
is dynamics usi ng the smoothed particle hydro- 
(SPH) method <MonagharJll992h . The envelope 



1 These parameters were chosen to roughly represent the inferred 
properties during the late spiral-in phase for merger models of 



SN 1987A, where the core fraction represents the immersed binary 
core, consisting of the core of the primary and the spiraling-in 
companion. 
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density is sampled with 10 5 particles using a Monte-Carlo 
method followed by isentropic relaxation to reduce numeri- 
cal noise llLucvl 19771) . The code units are M, the mass of the 
envelope and R, the initial (non-rotating) stellar radius (the 
total stellar mass including the core is 5M/3). This implies 
that time is measured in the code in units of y 1 R? /GM and 
velocity in units of yj GM/R. Note that this allows simple 
re-scaling of the results presented in this paper. 

To parametrize the spin-up of the envelope and the 
energy deposited by the spiral-in, we define two parame- 
ters a and j3, where a = AE/Eb is the ratio of the en- 
ergy deposited to the binding energy of the envelope and 
f3 = L/y GM 3 R is a dimensionless measure of of the enve- 
lope angular momentum following the early in-spiral of the 
secondary. 

To spin up the envelope we add angular momentum 
on a dynamical timescale using the following recipe: during 
every fixed timestep At = 0.025, the rotational velocity of 
each particle is incremented by an amount 



Ay, = riAtt, provided v% < y ~4>i, (6) 



i.e. as long as the velocity remains sub-Keplerian. The angu- 
lar velocity increment is Af2 = 0.0093 of the critical velocity 
at the surface of the non-rotating envelope (r = 1 in code 
units). If at any time the particle velocity Vi reaches the lo- 
cal Keplerian velocity, we set Avi — thereafter. This leads 
to solid body rotation in the inner envelope and a slightly 
rising specific angular momentum profile in the outer enve- 
lope (see Fig. 0. The spin- up phase is terminated when f3 
reaches 0.235, 0.588, 0.817, respectively, for the three cases 
we consider in this paper (note that no particles become 
super-critical; see Figure and Table 0. These three val- 
ues correspond to envelope angular momenta of 2.3, 5.7 and 
8.0 x 10 54 gem 2 s _1 which is comparable to the available or- 
bital angular momentum (equation 0. 

To simulate the deposition of the energy and the rapid 
heating of the envelope during the plunge-in phase, we then 
add entropy to the i nner envelope \r < 2/15R , from cal- 
culations discussed i n llvanova fc Podsiadlowskil teOOSt) and 
iPodsiadlowskil J200ll) ] . Initially, we consider the case where 
the energy is deposited immediately after the envelope has 
been spun up. The response of the envelope is then followed 
for 10 - 15 dynamical timescales after the instantaneous en- 
ergy deposition, at which point all the ejected particles are 
to good approximation on ballistic trajectories. 

Since we do not follow the evolution of the spiraling-in 
binary components, we do not encounter resolution prob- 
lems when the orbital sep aration becomes com parable to 
the SPH smoothing length IILivio fc Sokerl ll988 % ): our model 
is mainly limited by physical approximations (such as the 
equation of state; the lack of energy transport) rather than 
the numerical resolution (see Appendix A) - except near the 
surface. The steep density profile at the surface is poorly re- 
solved by SPH particles, and a well known problem of SPH 
models in this context (as, e.g., seen in supernova models) 
is an over-estimate of the mass contained in low-velocity 
material. 



Table 1. Properties of the three rotating envelopes immediately 
before the energy deposition for the zero-delay case. All values in 
cgs units are for M env = 12 Mq, R = 1500 Rq appropriate to the 
merger model for SN 1987A. 



T/W 


0.039 


0.117 


0.151 


13 = L/VGM :i R 


0.235 


0.588 


0.817 


Angular momentum 10 54 gcm 2 s _1 


2.3 


5.7 


8.0 


Mean angular velocity 10 -8 s — 1 


3.7 


2.2 


1.1 


Flattening ,Req/i?polar 


1.7 


4.3 


6.8 


Rotation velocity (km/s) 


42 


21 


15 


Fraction of critical rotation 


0.96 


0.69 


0.61 


Binding energy (10 47 ergs) 


-6.2 


-5.4 


-4.5 



Table 2. Properties of the rotating envelopes immediately be- 
fore the energy deposition, for L = 0.588\/ GM 3 R (top) and 
L = 0.817V GM S R (below) and for three values of the time delay 
following the spin-up of the envelope. 



<dciay (years) 


0.67 


1.35 


2.69 


T/W 


0.090 


0.078 


0.071 


Angular momentum 10 54 gcm 2 s _1 


5.7 


5.7 


5.7 


Mean angular velocity 10 -8 s — 1 


1.46 


1.07 


0.71 


Flattening ,Req/i?polar 


5.2 


5.4 


7.0 


Rotation velocity (km/s) 


18 


16 


13 


Fraction of critical rotation 


0.67 


0.63 


0.57 


Binding energy (10 47 ergs) 


-5.2 


-5.1 


-5.0 


T/W 


0.121 


0.104 


0.088 


Angular momentum 10 54 gcm 2 s _1 


8.0 


8.0 


8.0 


Mean angular velocity 10 -8 s — 1 


0.79 


0.61 


0.42 


Flattening i? e q/i? P olar 


7.5 


8.3 


9.6 


Rotation velocity (km/s) 


14 


13 


11 


Fraction of critical rotation 


0.58 


0.56 


0.52 


Binding energy (10 47 ergs) 


-4.3 


-4.2 


-4.1 



3 RESULTS AND DISCUSSION 

As shown in Figure 1, the addition of substantial angular 
momentum strongly distorts the envelope (see, in particu- 
lar, the 8 x 10 54 gcm 2 s _1 calculation and Table 0. These 
profiles are similar to the (n, n') — (1.5, 0) sequence of invis- 
cid po lytropic models calculated in lBodenheimer fc Ostrikerl 
( 1973), although our envelopes are more extended for a given 
value of T/W (the ratio of kinetic to potential energy) due 
to some viscous heating. The outer envelopes contain 1.7 
and 4M for (3 = 0.588 and 0.817, respectively, and both 
have a temperature of ~ 10 4 K. 

Following the energy deposition, matter is ejected in a 
very anisotropic way depending on the rotational distortion 
and the amount of energy deposited (i.e. depend on a and /3; 
see Fig. [3J . Generally envelope material is first ejected (i.e. 
reaches escape speed) in the polar direction. At low rotation 
(/3 = 0.235, not illustrated) mass is also ejected in the equa- 
torial plane at early times, and the distribution is more or 
less spherically symmetric. For larger values of /3, envelope 
ejection in the equatorial plane is suppressed for low values 
of a or enhanced for large values of a. In particular, no mass 
is lost in the equatorial plane if the energy fraction a is less 
than some critical value a c which increases with increasing 
/3. If no matter is ejected in the equatorial plane, we find a 
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a=0. 25,0 = 0. 817 a = 0.33,/S = 0.817 a=0.40,/S=0.817 a = 0.50,/S = 0.817 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

COS0 COS0 COS& COS0 



Figure 3. The final amount of mass ejected (left scale) as a function of cos0 where 9 is the polar angle for different values of the energy 
deposition and different values of the angular momentum of the envelope (the values of a and fi are given above each panel; a = 1 
corresponds to an energy of 5.4 X 10 47 ergs (for /3 = 0.588) and 4.5 X 10 47 ergs (for /3 = 0.817) while /3 = 1 corresponds to an angular 
momentum L cnv = 9.7 X 10 54 gem 2 s _1 . The central solid curves show the median velocity of the material (right scale) ejected as a 
function of polar angle, and the upper and lower curves give the range of velocities that includes 50 % of the material. 
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Figure 4. The (/3, a) plane at t = 10, for the zero delay case. The line separates models with (above) and without (below) equatorial 
mass ejection. Labelled models correspond to models 1 and 2 discussed in the text. 



strong mass excess at mid-latitudes due to shock focusing 
by the extended envelope (as discussed further below). 

In Sections 13 . 1 1 and 13 . 21 we discuss two typical calcula- 
tions with a = 0.33, f3 = 0.82 and a = 0.35, (3 = 0.66, 
respectively. In the j3 = 0.82 case, the critical energy frac- 
tion is a c — 0.39 (from Figure and therefore no material 
is ejected in the equatorial plane. In both these calculations 
energy is deposited immediately after the spin-up phase. 

In Section 13.31 we discuss the mass loss geometry due 
to a variation of a and /3 for different time delays between 
the spin-up and energy deposition phase (no time delay and 
a delay of 1.35 yr, respectively; see also Table 2). The dy- 
namical time delay is a free parameter in our model which 
could be constrained by more detailed models that follow 
the spiral-in explicitly. 



3.1 Model 1: a = 0.33, (5 = 0.817 

In this model, the envelope is significantly distorted by the 
rotation (column 3 of Table 0, and most mass is ejected 
at a latitude between 30 and 40°. Despite the reduction in 
the effective gravity close to the equator, very little mass is 
ejected there. The absence of ejected material in the equa- 
torial plane is relatively easy to understand from the strong 
rotational deformation, since the shock wave cannot pene- 
trate the outer envelope in the equatorial direction. 

However, it is not quite so obvious how the mass en- 
hancement at mid-latitudes arises. To illustrate the origin 
of this enhancement, we plot in Figure |5] the evolution of 
the angular momentum in the z-direction relative to the to- 
tal angular momentum for selected particles, which eventu- 
ally escape, as a function of distance from the rotation axis 
during the early ejection phase. The ratio L z j L parametrises 
the inclination of the orbit of a particle, and its change shows 
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Modell 




sqrt(x 2 +y 2 ) 

Figure 5. L z /L as a function of distance from the rotation axis 
for selected (ejected) particles in model 1 to illustrate the devia- 
tions from ballistic motion, where L z is the angular momentum in 
the z direction and L the total angular momentum of the particle. 
Stars and circles show particle positions at t = 2 and t = 15 (in 
code units), respectively. The vertical axis shows L z /L; equatorial 
orbits are close to the top of the figure. 



Modell 




1 2 4 6 8 10 



Figure 6. Particle snapshot at t = 3 showing the deflection of 
the ejected particles (circled) from the outer envelope. 



the effects of the strong pressure gradients that act during 
this phase. Particles initially ejected at mid-latitudes move 
poleward shortly after ejection, reach a minimum in L z /L 
before t — 3 (e.g. the blue trajectory in Figure Qj] which 
reaches a minimum at L z /L = 0.5) and then evolve towards 
a ballistic orbit with a lower inclination. 

During the strong poleward deflection at t = 2, a bow 
shock forms ahead of the massive outer envelope (see Fig- 
ure EJ, which deflects particles away from the equatorial 
plane. Hence the mass loss poleward of = 30° changes 
from isotropic at t — 2 to strongly peaked at t = 3, where 8 
is the co-latitude. 
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Model2 




5 10 

sqrt(x 2 +y 2 ) 

Figure 7. Similar to Figure 2 for model 2. Stars and circles show 
particle positions at t = 2 and t = 15 (in code units), respectively. 
The vertical axis is L z /L. 



The second deviation in Figure can be understood 
by the reflection of particles from slow-moving material at 
higher latitudes, which was ejected slightly earlier. This 
leads to a compression of the outflow at mid-latitudes 
{L z /L ~ 0.65), increasing the density enhancement. After 
t = 5 the particles follow nearly radial trajectories, although 
the trajectories are not yet quite ballistic. By t = 10 the flow 
is almost completely ballistic (as indicated by the fact that 
L z /L remains constant thereafter). 



3.2 Model 2: a = 0.35, (3 = 0.66 

The principal effect of the reduced angular momentum of 
the envelope is that the shock wave eventually reaches the 
surface in the equatorial region, and that some material is 
ejected there at low velocities (similar to the a = 0.33 and 
13 = 0.588 model in Fig. 3 but with a larger excess at mid- 
latitudes). The total ejected mass at t = 10 is 0.51 Mq, of 
which 0.04 is ejected per unit solid angle in the equa- 
torial plane, a value that is close to the value one would 
expect if the ejected mass had spherical symmetry. Initially 
the flow is compressed due to th e Bjorkman-Cassinelli ef- 
fect toiorkman fc Cassinellilll 993), but a strong density en- 
hancement of one to two orders of magnitude does not form 
due to the fact that our envelope is extended and mainly 
supported by thermal pressure. 

The mid-latitude enhancement contains a similar 
amount of mass as Model 1 discussed above, although its 
latitude is closer to the equator. As can be seen in Figure |7| 
the flow is qualitatively similar to Model 1 (Figure 0, but 
the shock interactions occur at lower latitudes since the star 
is less rotationally flattened. The velocity profile is similar 
in the polar region, since almost the same amount of energy 
was deposited in the envelope. 

As in Model 1, the initial polar enhancement is caused 
by the absence of ejected matter elsewhere at early times. 
This leads to a large asphericity factor as a result of the small 
solid angle into which mass is ejected. Once mass has been 
ejected at lower latitudes, the 'enhancement' disappears. 
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= 0.25,(3 = 0.817 



a = 0.33,/S = 0.817 



= 0.40,^=0.817 



= 0.50,(3 = 




0.2 0.4 0.6 O.f 

cosf? 




0.2 0.4 0.6 O.f 
COS0 




2 0.4 0.6 Of 
COS0 




0.2 0.4 0.6 O.S 
cosf? 



Figure 8. Properties of the ejected matter, similar to Figure [3] for the case where the energy is deposited 1.35 years (1.6 dynamical 
timescales) after the spin-up of the envelope. 



Figure 9. Total ejected mass for the calculations listed in Table |3] after ~ 8.2 years. Darker shading corresponds to longer delay times. 



/3= 0.588 




0.3 0.4 0.5 

Relative energy deposition a 



/?= 0.817 




13 0.4 0.5 

Relative energy deposition a 



3.3 Parameter study 

In general, the mass ejected depends on three principal pa- 
rameters, the amount of energy deposited (AE), the total 
angular momentum in the envelope (L) and the time delay 
between the spin-up of the envelope and the deposition of 
the energy (tdeiay), which depends on the timescale of the 
initial spiral-in phase. Variation of these parameters leads 
to changes in i) the total ejected mass, ii) the presence or 
absence of ejected material in the equatorial plane, and iii) 
the strength of the enhancement at mid-latitudes. The peak 
mass flux also moves to slightly higher latitudes with in- 
creasing rotation. 

In Figure we show the results of a systematic param- 
eter study. Each panel shows the geometry of the ejected 
mass as a function of cos (9, where 9 is the polar angle, once 
the ejecta are expanding ballistically on radial trajectories. 
The histograms show the mass ejected divided into bins of 
constant dcosf (each bin subtends 7r/10 sr), and the indi- 
vidual curves show the velocity distribution at each angle 
(the central curve gives the median velocity, and the upper 
and lower curves give the velocity range which includes 50 % 
of the ejected matter). The velocity and angular momentum 
scale according to 



n C r ~ 39kms 



x / M 1500_R Q x 1 ~ 
\ 12M Q ~R~ 



(7) 



^^arxio-gcm^s-^) 372 ^) .(8) 

The peak at mid-latitudes, which is only present for 
sufficiently distorted envelopes with j3 > 0.5, remains up 
to a = 1. The mass excess at mid-latitudes becomes more 
pronounced as the angular momentum increases. The ve- 
locity of the equatorial material, when present, is typically 
a factor of 3-5 lower than the velocities at mid-latitudes, 
though this difference is reduced as the deposited energy is 
increased. The critical a below which no mass is lost in the 
equatorial plane increases with rotation rate (see Figure HJ, 
since the envelope is more extended. 

The trends remain very similar, when we introduce a 
time delay between the spin-up phase and the energy depo- 
sition phase (see Figure|HJ , except that the loss of material in 
the equatorial plane is now further impeded by the massive 
extended envelope, and the mid-latitude enhancement may 
be much stronger (see, for example, the a = 0.5, P = 0.817, 
tdciay = 1-35 yr case). The peak has moved to slightly higher 
latitudes compared to the zero delay case, but this is a small 



1/2 
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Table 3. Total mass ejected at t = 10t t j Lyn ~ 8.2 years after the entropy deposition in the inner envelope, in solar masses. T is the time 
delay in years following the spin-up phase. 



a = 


E 


8 - . L , - 0.235 
T = 


T = 


3 = 
T = 0.67 


0.588 
T = 1.35 


T = 2.69 


T = 


3 = 
T = 0.67 


0.817 
T = 1.35 


T = 2.69 


0.25 




0.24 


0.17 


0.14 


0.17 


0.21 


0.18 


0.15 


0.16 


0.18 


0.33 




0.44 


0.55 


0.32 


0.36 


0.42 


0.34 


0.32 


0.34 


0.38 


0.4 




0.83 


1.12 


0.51 


0.55 


0.63 


0.68 


0.49 


0.51 


0.58 


0.5 




1.33 


2.15 


1.10 


0.91 


1.00 


1.67 


0.79 


0.81 


0.89 



effect. During the time delay the outer envelope expands by 
a factor ~ 1.5 with a corresponding increase in the critical 
energy for equatorial mass ejection. Hence increasing the 
energy drives more mass from each surface element but has 
little effect on the geometry. The total mass ejected for each 
value of the energy deposition a and time delay is listed in 
Table 3 and illustrated in Figure 9. 

3.4 Summary 

Mass ejection during a common-envelope phase leading to 
the complete merger of a binary system preferentially occurs 
at mid-latitudes due to shock reflection from an outer enve- 
lope containing a significant fraction of the initial orbital an- 
gular momentum. If the energy deposited is less than about 
one-third of the binding energy of the envelope no material 
is lost in the equatorial plane, in contrast to earlier models 
which co nsidered the higher energy case appropriate for CE 
ejection (L ivio fc Sokerlll988l : ISandauist et al.lll998T) . 

Alt hough our results can not be directly compared to 
those o flSandauist et all Jl998f) since we do not model the in- 
spiral of the secondary, a number of similarities are apparent. 
The peak energy dissipation occurs when the orbital sepa- 
ration is approximately one-tenth of its initial value and the 
red giant envelope is already aspherical which is consistent 
with our prescription. Most of the mass loss (~ 0.3Af cnv ) is 
concentrated in the equatorial plane in this more energetic 
case. The equivalent energy parameter is a = AE/Eb ~ 1.2 
whereas we consider the range a = 0.25 — 0.5. 



4 APPLICATIONS TO SN1987A AND 
SHERIDAN 25 

Many non-spherical nebulae are axially symmetric, which 
has been interpreted as evidence for rapid rotation, pos- 
sibly as a result of binary interactions. A notable exam- 
ple is the mysterious nebula surrounding the supernova 
SN 1987A, which consists of three roughly parallel rings, 
one centered on the supernova and t he other two displace d 
by ~ 1 arc-second to either side feurrows et al.l Il995l) . 
The supern ova itself was anom alous in several other re- 
spects (see iPodsiadlowskil Il992j and references therein), 
which are most consistent with a binary merger some 
20,000 years before the supernova event jPodsiadlowski 
ll992l:|Podsiadlowski fc Ivanoval2003T) . In particular, the blue 
supergiant progenitor and the chemical anomalies in the 
inner ring can easily be explained as the result of the 
dredge-up of core material in t he final stage of the merger 
lllvanova fc Pod siadlowski 200j|). 



The fast wind (M ~ 10~ 7 Af Q yr" 1 , ~ 500km/s) 
of such a blue supergiant will sweep up and enhance any 
structures already present in the ejecta, leading to an axially 
symmetric but highly aspherical nebula. 

4.1 Late case B/C merger: SN 1987A 

During late He shell burning the primary will be a red su- 
pergiant with a critical surface rotation velocity of 

-if M 1500^0 \ V2 
^~ 45kmS [ l5M Q R ) (9) 

which is comparable to the velocities observed in the 
SN1987A nebula (see Table gj. The latitude-dependence of 
the ejected material of Model 1, shown in Figure [TTJ1 is char- 
acterised by a strong enhancement at mid-latitudes, while 
no material is lost in the equatorial region. The following 
features of the nebula may therefore be understood: 

(i) Strong mass enhancement in the outer rings. The 
outer rings (ORs) are a real density enhancement (lOOx 
the ambient value) and are not simply due to limb b righten- 
ing of an hourglass structure feurrows et aljFl995h . In our 
model the ORs result from wind-driven pressure gradients in 
the seed structures which directly result form the anisotropic 
ejection of material during the merger phase. Previous mod- 
els based only on equatorial density enhancements in the 
pre-existing material ha ve been unable to expl ain the high 
density in the ORs (e.g. iMartin fc Arnettj|l995l) . 

(ii) Displacement of the outer rings relative to the inner 
ring. We favour a model in which the inner ring originates a 
few 1000 years after the merger event, in a rotation-enforced 
outflow during contraction on the post-merger blue loop 
(Heger & Langer 1998; also see Collins et al. 1999). The rel- 
ative displacement of the outer rings can be understood if 
the mass ejection during the merger event itself was slightly 
asymmetric, perhaps due to a non-axisymmetric pulsational 
instability in the envelope, which gives the ejecta a veloc- 
ity of ~ 2kms _1 relative to the merger remnant. Hence the 
wind-driven pressure gradients are no longer axisymmetric 
and the planes 2 in which the outer rings lie will be slightly 
inclined with respect to the plane of the inner ring. This 
would explain both the offset of the outer rings and their 
shape, which is noticeably non-elliptical in projection. 

(iii) North/South asymmetry. Since the planes of the 
outer rings are no longer parallel to one another, the South- 
ern outer ring is observed closer to face-on in projection than 
the Northern one. 

2 The outer rings are still approximately planar in this case. 
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Table 4. Observed properties of circumstellar material around massive stars. Bracketed periods denote systems 
which are too wide to undergo binary interactions. The properties of the material in the equatorial torus are given 
in the last four columns, where R, v and t give the radius, expansion velocity and dynamical age, respectively. 



Star 


Spectral type 


Period 


Luminosity 
L/Lq 


Equatorial torus 
Density (cm -3 ) or mass R (10 17 cm) 


V (km/s) 


t (yrs) 


Sk -69° 202 


B3 la 




10 5.2 


n e ~ 2 - 5 X 10 4 


6.23±0.08 a 


10.3 


19,000 


Sher 25 


B1.5 lab 




1Q 5.9 


0.01 - 0.1M© 


6.2 


20 


6,600 


RY Scuti 


06.5 + 09.5 6 


11.3d 


~ 10 6 


1M S ? 


0.2 


50 


130 


r\ Car 


? Ia 


(5.5yr) 


10 6.7 


3 - 15M^ 


0.7 


50 


~ 1,000 


R4 d 


B[e] + A 


(21.3yr) 


10 5 + io 4 - 2 


? 


? 


? 


~ 10, 000 



Note: Data taken fro m HPanagial <2004l) . ISmith et all fcOOllK iMorris et all ll99Sft : ISmith et all i2003al) and 
' jZickgraf et al 1 <199rt . 





Figure 10. Mass enhancements in the ejecta flow, corresponding to merger models for SN 1987A (left-hand panel, Model 1 from 
section 3.1 with a = 0.33 and /3 = 0.817) and Sheridan 25 (right-hand panel; Model 2 from section 3.2 with a = 0.35 and f3 = 0.665). 
The solid curves give the locations that contain 50% of the mass ejected at a particular solid angle at the time as indicated (in code 
units). 



4.2 Case B merger: Sher 25 

Mass loss during common-envelope evolution may explain 
the broadly similar structures seen around other luminous 
stars, listed in Table|l] Of these, the nebula around the B1.5 
supergiant Sheridan 25 shows the most compelling similari- 
ties, since it has an equa torial ring of at least 0.01 — O.IOMq 
jBrandner et alJ Il997bl) and polar lobes, each containing 
~ 0.25Mq. It is in a post main sequence, though proba- 
bly pre-red s upergiant, evolutio nary state (with N/C ~ 26, 
N/O = 0.36. ISmartt et alJliool) . The dynamical age of the 
nebula has been estimated to be around 6,000 years. The 
observed velocities suggest an envelope radius of ~ 300-Rq 
at the time of ejection, corresponding to a merger during 
the Hertzsprung gap crossing. 

Although the density structure of the common enve- 
lope will differ from that of a 7 = 5/3 condensed polytrope, 
the results may still be applicable if the envelope is aspheri- 
cal. We suggest a model with equatorial mass loss during 
the merger, such as Model 2 discussed above (see right- 
hand panel of Figure HUB , since the envelope cannot store 
enough angular momentum to generate a significant post- 
merger equatorial outflow. Hence an asymmetry during the 
merger will displace both the equatorial and polar material 
from the site of the merged st ar. The equatorial ring is offset 
from Sher 25 by 0.05-0.1 pc ferandner et alJll997al) . which 
is consistent with this model. 



4.3 Conclusions 

The three principal anomalous features of the supernova 
SN1987A, viz. its blue supergiant progenitor, its over- 
abundance of certain elements, notably He, and the presence 
of highly structured circumstellar material, are all consistent 
with a binary merger some 20,000 years before the explosion. 
In this paper, we have demonstrated how density enhance- 
ments at mid-latitudes arise during mass ejection from a 
rotationally distorted star. Subsequent interaction with the 
fast wind of the blue super giant prior to the supernova (cf. 
iBlondin fc Lundavistl ll993) then leads to the formation of 
the outer rings wi th a density enhancement of a fa ctor of 150, 
in calculations of lMorris fc Podsiadlowskil j200fj) which will 
be further discussed in a future paper. 

Similarly, the nebula around Sher 25 may be explained 
by a binary merger following a CE phase during the cross- 
ing of the Hertzsprung gap by the primary. One notable 
difference in this case is that the equatorial ring likely orig- 
inates during the merger, which is consistent with observa- 
tions showing that the centre of the equatorial ring is dis- 
placed by some 0.05 pc from Sher 25. Future observations 
of the rotation rate of Sher 25 would help to confirm this 
model. 

Although the nebula around the more massive sys- 
tem RY Scuti (09.5 + 06.5) appears similar, its evolution 
is probably somewhat different. Data from the Keck tele- 
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scope and the HS T, discussed in ISmith et alJ il999T) and 
ISmith et all J200ll) . show a massive equatorial dust torus 
and two narrow rings symmetrical about the equatorial 
plane, with a dynamical age of 120 years. The massive torus 
probably originated during thermal timescale mass outflow 
from the outer Lagrangian L2 point which is still occurring 
today, albeit at the much lower rate of ~ 5 x 10~ 5 Mq yr _1 . 
Subsequent deflection of the fast wind i n a manner s ome- 
what analogous to the horseshoe model o f lSokerl dl999f> may 
expl ain the origin of the two narrow rings. 

IPasauali et al.l i2000t) have suggested that the B [e] com- 
ponent of the spectroscopic binary R4 in the LMC could 
be the result of a Case B/C merger, i.e. the system would 
originally have been a triple system, where the companion 
A star now serves as an astronomical clock and indicates 
that the present primary has lost > 40% of its z ero-age 
main-sequence (ZAMS) mass dZickeraf et alJll996l) . Since 
the companion is already a post-main sequence A super- 
giant, the pre-merger components cannot be much more 
mass ive, e.g. 12 — 1 5 Mr for the original primary. In addi- 
tion, IPasauali et al.l (l200CTt find evidence of CNO-processed 
circumstellar material with a dynamical age of ~ 10 4 years 
which is consistent with ejection during a binary merger. 
Direct images of the nebula would be of great benefit to 
understanding this system. 

Bipolarity is also common in observations of planetary 
nebulae. Asymmetric mass loss during a common-envelope 
phase (with or without a merger) provides physical motiva- 
tion for the equatorial density enhanc emen t functions pro- 
posed bv llcke. Ba.lick fe Prestonl Jl999h and lLuo fe McCravl 
il99lfl (see also lFranklll999l) We speculate that the ho- 
munculus nebula around r\ Carinae may also have originated 
during a co mmon-envelope phas e since its kinetic energy is 
~ 10 50 ergs jSmith et alJl2003al) . comparable to the lumi- 
nous energy of the 1840 — 1860 outburst. Both the mass 
loss rate in the stellar wind (1.6 x lO _3 M0yr _1 ) and its 
latit ude dependence suggest rapid rotation of the central 
star (ISmith et alJl2003bt Ivan Boekel et alj20ollAerts et alJ 
12004^1 . In the merger scenario r] Carinae was originally a 
triple system in which the closer components (P ~ 30 d) 
merged 150 yr ago to leave the present companion in an ec- 
centric 5.5-yr orbit. 
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COS0 COS0 

Figure Al. Comparison between two calculations of Model 1, 
with 2 X 10 5 particles (left-hand panel) and 4 X 10 5 particles 
(right-hand panel). 



This paper has been typeset from a TfjjX/ IflTjjX file prepared 
by the author. 



APPENDIX A: RESOLUTION STUDY 

Figure lATl shows that increasing the number of particles 
by a factor of two does not significantly change the geometry 
of the ejecta. This implies that our calculations have con- 
verged numerically. The calculation is limited by the phys- 
ical approximations such as the assumption of a polytropic 
equation of state. 
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